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Abstract. The quantum Zeno effect describes the inhibition of quantum evolution by 
frequent measurements. Here, we propose a scheme for entangling two given photons 
based on this effect. We consider a linear-optics set-up with an absorber medium whose 
two-photon absorption rate ^27 exceeds the one-photon loss rate 1^1^. In order to reach 
an error probability Pcrror, we need ^17/^27 < "^Pcttox I '^'^ ^ which is a factor of 64 better 
than previous approaches (e.g., by Franson et al ). 

Since typical media have ^27 < C17J '^6 discuss three mechanisms for enhancing 
two-photon absorption as compared to one-photon loss. The first mechanism again 
employs the quantum Zeno effect via self-interference effects when sending two photons 
repeatedly through the same absorber. The second mechanism is based on coherent 
excitations of many atoms and exploits the fact that ^27 scales with the number of 
excitations but ^i-^ does not. The third mechanism envisages three- level systems where 
the middle level is meta-stable (A-system). In this case, fi-y is more strongly reduced 
than ^27 and thus it should be possible to achieve £,2-y/£,i-t ^ 1- 

In conclusion, although our scheme poses challenges regarding the density of active 
atoms/molecules in the absorber medium, their coupling constants and the detuning, 
etc., we find that a two-photon gate with an error probability Pcrror below 25% might 
be feasible using present-day technology. 
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1. Introduction 

The idea that quantum mechanics or, more precisely, quantum systems could be used 
for (quantum) computation [IJ-E] , led to the advent of a whole new research area known 
as quantum information processing. Important theoretical progress has been made in 
this field, showing that quantum information processing exceeds classical information 
processing in a fundamental sense [2ll3H6]. 

Accordingly, strong efforts are undertaken to physically realize quantum computers, 
following such different approaches as ion traps, quantum dots, Josephson junctions, 
nuclear spins and quantum optics, all with their respective assets and drawbacks [7]. 
Picking up one of the promising routes to quantum information processing, this work 
is connected with a quantum optics approach. In this approach, the information is 
encoded in the polarization or spatial degree of freedom of optical modes, i.e. qubits 
are physically implemented as single photons. 

It is well known that photons possess many properties which make them very 
suitable candidates for quantum information processing. They can be well controlled 
and manipulated, as well as created and measured, even down to the single-photon level 
P4TT]. Furthermore, because of their long decoherence time, photons can propagate over 
relatively long distances without significantly coupling to the environment. Moreover, 
when thinking of the quantum circuit model, see e.g. [12], the necessary single-qubit 
quantum gates (such as Hadamard and phase shift gates) exist for photons and are 
easily implemented. If an adequate two-photon entangling gate (such as a CNOT-gate) 
with low enough error-probability existed, one would have a universal set of gates, from 
which every quantum algorithm could be constructed. Unfortunately, the weak coupling 
of photons to the environment is also related to their main drawback - it is very hard 
to make two photons interact, and thus hard to build a two-photon entangling gate, as 
the CNOT-gate. Typically, before two photons interact by means of some non-linear 
medium, at least one of them is absorbed. This motivates the idea to turn the problem 
around and to actually exploit the absorption in order to make photons interact. 

An indirect way of doing this is realized by the Knill-Laflamme-Milburn (KLM) 
proposal [13] , which initiated the field of linear optical quantum computation. In the 
KLM paper it was shown that scalable quantum computing is theoretically possible 
using only linear optical elements (like beam-splitters and phase shifters), single photon 
sources and detectors. In this approach, interactions are induced (probabilistically) 
by means of entangled ancilla photons on which measurements are performed. Linear 
optical quantum computation is being studied in great detail theoretically (e.g., [T41 - I29] ) 
and experimentally (e.g., [30]-ST]), and is currently one of the most elaborate schemes 
for optical quantum information processing. A major drawback is the vast amount of 
resources (ancilla photons, entanglement) which is - despite all theoretical progress - 
still needed when scaling up this scheme. 

A more direct way of exploiting the absorption of photons is to employ the quantum 
Zeno effect which describes the slow-down or even inhibition of quantum evolution by 
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repeated measurements [l^Hlg] . Imagine a quantum particle in a double- well potential 
initially confined to the right well, for example. After some time T, it would tunnel to 
the left well (and then back, etc). However, if we measure the position of the particle 
frequently, i.e., after very short time intervals At <^ T, it does not tunnel since each 
measurement projects the quantum state back to the right well. This is the basic 
picture of the quantum Zeno effect. Since the absorption of a photon in some medium is 
equivalent to measuring its position, a strong enough absorption probability can actually 
prevent the photon from tunnelling/propagating into this medium |49j . 

The idea to utilize the quantum Zeno effect to build a two-photon entangling gate 
was first brought up by Franson and co-workers [50] and further developed in [511 - 155] . 
In their set-up, the two photons (qubits) are guided within two fibre cores, which 
are coupled by their evanescent fields. When doping the fibre cores with two-photon 
absorbing atoms/molecules, unwanted two-photon amplitudes are suppressed due to 
the quantum Zeno effect. One advantage of this approach is that it requires much less 
resources when compared to linear optics quantum computation (KLM-style) schemes. 
However, an important issue of the Zeno gate is that while strong two-photon absorption 
is strictly required, one-photon loss leads to its failure. That is, the rate of two-photon 
absorption needs to be several orders of magnitude larger than the rate of one-photon 
loss for successful gate operation. 

In this paper, after a brief introduction on the quantum Zeno effect (section ll.ip . 
a new optical set-up for an entangling quantum Zeno gate is developed (section |2]). In 
contrast to earlier proposals for photon gates based on the quantum Zeno effect [501455] , 
the set-up presented in this paper is modified and so yields a significantly reduced error 
probability, see ( j28l) . which was derived considering finite two-photon absorption and 
non-negligible one-photon loss. 

Furthermore, a realization based on free-space propagation is proposed, i.e., without 
waveguides, resonators, or optical fibres - which may induce additional decoherence. 
Thus in section [3l an analysis of two-photon absorption and one-photon scattering is 
given for a single three-level atom in free-space. As our gate requires an absorbing 
material which features strong two-photon absorption compared to one-photon loss, 
we present three mechanisms (sections HI [5l and [T]) to enhance two-photon absorption 
compared to one-photon loss, which could be applied in free-space. Additionally, by 
inserting example values (sections |6] and [7]), we show that our scheme is not out of reach 
experimentally. We summarize and conclude in section [HI 

1.1. The quantum Zeno effect 

In our approach to entangle photons, we employ the quantum Zeno effect. The quantum 
Zeno effect describes the suppression of the time evolution of a quantum state by frequent 
measurements. It is based on a subtle application of the quantum theory of measurement 
processes [56]. Imagine a quantum particle in a double- well potential as depicted in 
Figure [TJ Let us assume that the quantum particle initially is confined in the right well. 
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Figure 1. Double- well potential with particle in the right well initially. 



1^(0)) = |R), (1) 

where a ket- vector |R) denotes the basis state where the particle is in the right well and 
|L) denotes the basis state where the particle occupies the left well. As the two wells of 
the double-well potential are coupled, according to the rules of quantum mechanics, the 
quantum particle starts tunnelling from the right well into the left well. After a certain 
time T, the quantum state will have evolved to 

l^(T)) = |L), (2) 

1. e., the particle will have tunneled entirely into the left well. Then the particle starts 
tunnelling back such that it is entirely in the right well at t = 2T etc. 

However, when dividing the tunnelling time T into ^ 1 equal-spaced time steps 
At = T/N, the quantum state after the first time step would be 

|^(At)) = cos(^)|R) + sin(^)|L), (3) 

and after the n-th timestep, the state would evolve into 

\ij (nAt)) = cos (n^) |R) + sin (n^) |L) . (4) 

So far, there were no measurements made on the system. Now imagine we measure 
the position of the particle frequently, i.e. after each time step At. This means 
that the quantum state ([3]) is projected back onto the initial state \ip (0)) = |R) with 
probability Pr (At) = cos^vr/ (2A^) or is projected onto the state |L) with probability 
Pl (At) = sin^ 71 / {2N). Furthermore, when measuring the position of the particle after 
each time step, the probability for the particle remaining in the right well is given by 

Ph, (T) = (At) = cos^^ (^) = 1 - + O i^-') ' (5) 
which approaches unity in the limit of large A^ ^ 1. With this probability the 
time evolution of the quantum state is completely supressed, i.e. the quantum particle 
stays in the right well instead of the left well, where it would have arrived without 
measurements. 

2. Two-photon entangling Zeno gate 

2.1. Two-branch gate 

In section II. the basic principle of the quantum Zeno effect was briefly discussed. 
In this section we will demonstrate how the quantum Zeno effect can be employed in 
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(b) 

Figure 2. Comparison of tlie dynamics of tlie particle with (b) and without (a) 
frequent measurements, (a) Without quantum Zeno effect, the quantum particle 
tunnels from the right well (t = 0) into the left well [t — T). (b) When measuring the 
position of the quantum particle after each time step At, the particle is still found in 
the right well for t ^ T. This is a simple example of the quantum Zeno effect. 

an otherwise linear optics set-up for building a two-photon entangling gate. For this 
purpose we will first study the underlying concept via a reduced version of our entangling 
gate, called the "Two-branch gate". On this basis we will be able to understand the 



functionality of the "Three-branch gate" (section 12.21) and to derive the dependency 




Figure 3. Sketch of the macro-structure of the two-branch gate. Horizontal solid lines 
(top and bottom) denote perfect mirrors and horizontal dashed lines indicate beam 
splitters. The target photon (slanting dark blue line) enters the upper branch. The 
control photon (slanting dark green dashed line) enters the lower branch (if present). 
Then both, the control photon and the part of the target photon which tunneled 
through the beam splitter into the lower branch, pass the two-photon absorbing 
medium (brown circle). Each segment (tilted light blue oval) consists of an absorbing 
medium and a beam splitter. 



The proposed quantum Zeno gate consists of perfect mirrors at the top and bottom, 
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weak beam splitters in the middle, and a two-photon absorbing medium in the lower 
branch, see Figure [31 The quantum state of the target photon is represented by a 2- 
vector ip = {x, y)'^ where x and y are the amplitudes for the target photon being in 
the upper and lower branch, respectively. Additionally, the control photon is either 
incident in the lower branch, or not, \0)fj. The polarization and/or frequency of 

the control photon is chosen such that it is perfectly reflected by the beam splitters and 
thus maintains its state in the lower branch for the whole gate. The target photon, 
in contrast, is affected by weak beam splitters with reflectivity cose and transmittivity 
sine, where e ^ 1. Their effect on the state ip can be represented by a 2x2 rotation 
matrix 



cos e sm e 
— sin e cos e 



(6) 



The absorbers in the lower branch reduce the amplitude by a factor e~'^, where ^ takes 
different values depending on whether the control photon is incident, ^27, or not, ^1^. 
Therefore the whole operation of one segment of the gate on the target photon state 
reads 



(7) 



cos e sin e 




1 




cos e sin e 


e~'^ sin e e~'^ cos e 




e-« 




— sin e cos e 



As our gate consists of A^ of such segments, the output state is given by 

-I N 

COS e sin e 



out 



sm e 



cose 



(8) 



To explain how the gate works, we first consider the ideal case where the absorber 
exhibits perfect two-photon absorption, ^27 = 00, and non-existent one-photon loss, 
^1^ = 0. If we assume an input state where the target photon enters the upper branch, 
ipin = (Ij 0)"^, and there is no control photon (thus ^ = ^17 = 0), then the small rotations 
(E]) of the target photon state simply add up 



(9) 



Physically, the target photon in the upper branch gradually tunnels into the lower branch 
(and would continue to tunnel back and forth afterwards) via the weak beam splitters. 
Choosing e to fulfill A^e = 7r/2, the target photon will have completely tunneled into 
the lower branch at the end of the gate, ifl'l^ = (0, —1)"^, as then A^ rotations by a small 
angle e add up to a 7r/2-rotation. 

However, if there is a control photon incident and thus ^ = ^27 = 00, the second 
row of the matrix ([71) is zero and therefore the output state ([H]) is given by 







TV 






cos A^e 


& = 


COS e sin e 




1 




— sin e cos e 









— sin A^e 



rout 





N 






N 

COS e 


COS e sin e 




1 


















(10) 



In this case, due to the quantum Zeno effect, the target photon stays in the upper 
branch. Physically, each beam splitter transmits a small fraction sin e of the amplitude 
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from the upper branch into the lower branch, which immediately gets canceled by the 
two-photon absorbing medium. Hence, the amplitude in the upper branch is reduced 
by a factor cose for each beam splitter. In the limit of large N, ijj^l^ 
approaches 1pl'l^ 
with O (AT-i) 



cos 



TV 



e,Of 

'1,0)^, where the error due to finite N (discretization error) scales 



cos 



N 



TC 



2N ~ ~8N 

For ideal absorbers (^27 = 00, = 0) and N ^ 00, the illustrated two-branch- 
gate would deterministically entangle photons, as the branch where the target photon 
ends up is totally determined by the presence or absence of the control photon. These 
requirements, however, are clearly unrealistic when it comes to real, physical absorbers. 
In the following, we therefore investigate how the gate operation is affected when 
considering finite two- photon absorption (^27 < 00) and non- negligible one-photon loss 
(^1^ > 0). For arbitrary parameters ^ and A^, ([8]) can be evaluated to 



+ O (iv-2) 



1 



out 



2^7 



{(5+a1 - P^a'l) /2 



a 



N 



a 



N 



) 



sme 



( 

[(5+a 



N 



OL"l — a^) sin e 

N 



- /2 



(12) 



with abbreviations 



e ^ + 1)^ cos^ e 



4e-« 



(e ^ + 1) cose ± r 



(13) 



a± = 

j3± = (e~^ — 1) cos e ± r . 

It is insightful to throw a glance at the output state f[T^ for the one-photon case, 
expanded for small ^ = <^ 1 in first order, where ipm = (1; 0)^. It reads 



rout 



O (e?7) . 



(14) 



cos A^e (1 - ^1^, (A^ - cot e tan A^e) /2) 
-sinAre(l -67 (Ar + l)/2) 

As was shown before (discussion below ([H])), the absolute value of the amplitude in the 
lower branch at the end of the gate is 1 in the ideal case, where = 0. Also in the 
general case, ,^1^ > 0, the absolute value of the second line in f lT^ needs to be to close 
to 1 in order to approximate the desired gate operation. This, however, is only possible 
when scales with 1/A^, such that A^^i-y <C 1. Moreover, ( !T^ shows that A^e = 7r/2 is 
still the optimal choice for e, even if the one-photon loss is non-zero. 



2.1.1. Error probabilities In both operating modes (with and without control photon), 
errors occur when considering non-ideal values for ^1^, ^27 and finite A^. We are going 
to derive the respective error probabilities and discuss the requirements on the absorber 
when the effective error probability of the gate should be below a certain error threshold 
p 

^ error- 

When the control photon is absent, the target photon could get lost (e.g. scattered) 
for ^ = > 0. This error corresponds to a situation where the target photon, which is 
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initially in the upper branch, ip^^ = (1, 0)-^, is not found in the lower branch at the end, 
thus 

pX = i-|(o,i)-v^;itf • (15) 

With control photon, finite N leads to a non-zero probability that the two photons 
actually get absorbed (discretization error, see discussion below ffTOl) ) and imperfect 
two-photon absorption, ^ = < 00, allows a non- vanishing amplitude in the ideally 
forbidden lower branch. In summary, the gate fails when the target photon is not found 
in the upper branch at the end 

2 



Pll. = 1 



(i,o)-v?;it • (16) 



Both error probabilities ( fTSi) . ( fT6l) are analytically well defined by ( [T2l) . but the exact 
expressions will not be quoted here. Instead we consider the limit N ^ 1, assuming 
that and ^27 scale with For this is necessary for proper gate operation, see 
discussion after dHj), and for ,^27 we need to choose the same dependence to obtain a 
quotient ^27/'Ci7 which is independent of (see below). Additionally, we demand that 
the one-photon loss rate is small even when multiplied with N, and that the two-photon 
absorption rate multiplied with is much larger than one (A^^27 ^ 1 ^ ^^17)- These 
limits are essential in order to obtain reasonable success probabilities. We then get for 
the error probability Pcj-ror i'^ '^^^^ without control photon (for the first non-vanishing 



order in 1/A^, see Appendix A.l ) 



PeLr = N^l, + O {N-') . (17) 

The error probability Parlor case with control photon reads 



error 
^2 



Usually both and ^27 scale linearly with the absorber length and therefore we can 
trade-off one error probability against the other. If we make the absorber longer, for 
example, both and ^27 increase by the same factor, which results in a lower P^ror 
but in a higher -Pgrror- 1^ define the effective error probability as the maximum 
of the error probabilities for both cases, the optimal absorber length is reached when 
-^error = -^crror' -Pgrror mouotouically increasing with the absorber length, while -P^ror 
is monotonically decreasing with the absorber length, see Figure HI Assuming that the 
quotient k, := ^27/^17 is fixed as a material property of the absorber, we arrive at optimal 
rates 

Reinserting this result into f[T7|) . f lTS]) yields the effective error probability for given k, 
optimized absorber length, and large A^ 

Perror = ^ + O {N-') . (20) 

\/2k 
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Figure 4. Plot of the error probabilities for the one-photon case (red) and for the two- 
photon case (blue) with k = 10^ and N = 10'^. The approximated error probabilities 
(IT71) . (dH]) are plotted as solid lines, while the dashed lines indicate the exact error 
probabilities p^ . ([TB| . ^27 — is proportional to the absorber length. The optimal 
absorber length is found at the intersection between the blue and the red graph. 



This relation can also be understood as a requirement on k in order to reach a certain 
error threshold Perror 

'5I7 error 

Before we examine ( I2T1) in detail, let us discuss the relation to earher proposals in [501I5T] . 
The proposed quantum Zeno gate in Figure [3] has substantial similarities to the quantum 
Zeno CSIGN-gate proposed by Leung and Ralph [51j based on the work of Franson et. 
al. [50]. However, the main difference here is that our set-up is not symmetric, i.e. 
the control photon always enters in the lower branch and stays there, and there are no 
absorbers in the upper branch. With the set-up described in this section, it is difficult 
to design a quantum CNOT-gate. 



2.2. Three-branch gate 

This obstacle can be overcome by making our set-up symmetric, i.e. adding a third 
branch at the bottom, such that the control photon always stays in the middle branch 
while the target photon can enter the gate in the top- or bottom-branch and tunnel to 
the opposite side (without control photon) or stay there (with control photon). Now the 
quantum state of the target photon is represented by a 3- vector ip = (x, y, z)'^, where x 
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Figure 5. Sketch of the macro-structure of the three-branch gate. Again, horizontal 
solid lines (top and bottom) denote perfect mirrors and horizontal dashed lines indicate 
beam splitters. The target photon (slanting dark blue line) enters the upper branch 
(initial state |0)) or the lower branch (initial state |1)). The control photon (slanting 
dark green dashed line) enters the middle branch (if present, i.e, if in the state 
Then both, the control photon and the part of the target photon which tunneled 
through the beam splitter into the middle branch, pass the two-photon absorbing 
medium (brown circle). If the target photon is not absorbed, it may continue to tunnel 
to the lower branch through the second beam splitter. Now, each segment (tilted light 
blue oval) consists of an absorbing medium and two beam splitters. 



stands for the amplitude in the upper branch, y for the amphtude in the middle branch 
and z for the amplitude in the lower branch. That is ipi^i = (1, 0, 0)"^ or ipi^ = (0, 0, 1)'^. 
The second line of beamsplitters is identical to the first and thus the operation of one 
segment of the gate on the target photon state reads 





cos e 



cos e sm e 
e~^ sin^ e 



— sm e 
e~^ cos^ e 
~^ cos e sin e 



1 
cos e — sin e 
sin e cos e 



1 









— sm e 
cose 



^ 
1 



cos e 
sine 




- sme 
cos e 






1 



Therefore the output state of the three-branch gate is given by 



out 



cose 



^ cos e sm e 
e~^ sin^ e 



— sme 
e~'^ cos^ e 



cos e sm e 





— sine 
cose 



N 



(22) 



(23) 



The relation between e and in order to have the target photon completely tunneled 
into the opposite branch at the end of the gate slightly changes into e = tt/ (-\/2A^) due 
to the additional beam splitter barrier. 



2.2.1. Error probabilities Now, without control photon (^ = ^i^), the error probability 
is given by the probability to not find the target photon, which is initially in the upper 
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branch, ipi^ = (1, 0, 0) , in the lower branch at the end of the gate 



p}:' = 1 



(0,0,1) -^1:1, . (24) 



2 



With control photon = ^27), as for the two-branch gate, the operation fails when the 
target photon is not in the upper branch at the end 

2 



p27 ^ I 



1,0,0) -iP^t ■ (25) 



Note that in both cases the error probabilities would be the same if the target photon 
alternatively starts in the lower branch and error probability definitions flMj) . (1251) change 
accordingly (setup is symmetric). 

V'out in (I23l) was calculated using Mathematica, but we will not quote the exact 
result here due to excessive length. Instead, we give the approximated results for the 
error probabilities flM|) . f l25|) . which were derived using the same assumptions as for the 
two-branch gate {N 1 and A^^27 ^ 1 ^ -^^17)- Interestingly, the results differ only 
slightly from the results described in section 12.1.11 (again, for first non- vanishing order 
in l/N, see Appendix A. 2 ) 

In analogy to section I2.1.H we are looking for the optimal two-photon absorption (^27) 
and one-photon loss (,^1^) rates, given a certain quotient k := ^7/^17. We arrive at 

67 = ' ^27 = V'^^ • (27) 



When reinserting fl27|) in fl26|) . we find that we achieve the same effective error probability 
f l20|) and the same constraint on the quotient between two-photon absorption and one- 
photon loss as for the two-branch gate, 

SI7 error 

The above requirement fl28l) represents a major experimental challenge since two-photon 
processes are typically much weaker than one-photon effects. There have been proposals 
to induce strong two-photon absorption in resonators, cavities, or fibres [57H60] . Of 
course, these ideas could be applied to our set-up in Figure [5] as well. But in the 
following sections, we shall focus on free-space propagation which offers some advantages 
as compared to wave-guides (but also has drawbacks, of course), and we are going to 
discuss the feasibility of the requirement ( l28|) on a concrete example. For an error 
threshold of Perror = 0.1, i.e. error probability below 10%, our gate would need a two- 
photon absorbing medium which features a ratio of k ^ 500. 



2.2.2. Comparison with Franson-Gate Let us briefly compare the error probabilities 
of our gate with the error probabilities of the set-up discussed in [SniEI]. When using 
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the same definitions for tlie error probabilities and applying the same limits, we obtain 
the following expressions for the Franson et. al. gate 

PX = 2iV67, (29) 
27r2 

Perlor = ^1, + + ^ (^"') " (30) 

It can be seen that the error probability for the Franson gate is higher for each case 
individually. Furthermore, the error probability for the two-photon case is always higher 
than the error probability for the one-photon case, which means that the optimal 
absorber length is given by minimizing f l30|) alone. Doing this (as always replacing 
^ = ^27/^17), we arrive at 

^^^ = 7^71^' ^^^ = ^71^' ^^^^ 

and thus we get an effective overall error probability of 

Looking at (152]) and (1281) . it is clear that the required value for k is a factor of 64 times 
smaller in our proposal compared to the Franson et. al. gate. This is mainly due to the 
fact that the absorbing medium is only in the middle branch in our set-up. 

2.2.3. Control photon loss Note that one-photon loss of the control photon was not 
taken into account so far. For a simple estimate of the additional error, assume as the 
control photon loss rate which should be small, A^^c 1- The error probability due to 
this effect is roughly 1 — (e"'*'') ~ 2A^^c- Adding this to our existing error probability 

Perror =^ + 2Ni, + {N'^) . (33) 

V 2/t 



20]), we get 



Assuming the worst case that the control photon loss rate is as high as the one-photon 
loss rate of the target photon, = ^17, we arrive at 

Perror = ^ (1 + 4) + O [N'^) , (34) 

and conclude that the error probability scales up by a factor of 5, and therefore the 
necessary ratio k in f l28]) increases by a factor of 25 (which is still less than the factor 
64 then Franson et. al. gate has in any case). For useful gate operation it is thus a 
necessity to make the one-photon loss of the control photon as small as possible. We 
will show later that it is justified to neglect the one-photon loss of the control photon, 
as we can make it much smaller than the one-photon loss of the target photon, ^ '^17- 
Summing up, the three-branch gate features the same probability of success as 
the two-branch gate, while offering the possibility to serve as a full-fledged quantum 
CNOT-gate. 
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3. Single three-level atom 




Both photons together are in 
El = uji + UJ2, but one photon 



As outlined in the previous section [21 the proposed quantum CNOT-gate requires a 
two-photon absorption rate which is several orders of magnitudes larger than the one- 
photon loss rate. In order to investigate if this requirement ( l28ll is satisfiable, we will now 
derive the two-photon absorption probability and one-photon scattering probability in 
the case of a simple three-level atom. The level scheme of the three-level atom is shown 
in Figure El it consists of three electron-eigenstates Is, 2p and 3s with their respective 
energy levels Ei, E2 and E^. An alternative level scheme is investigated in section [71 



Figure 6. Sketch (not to scale) of the level scheme, 
resonance with the transition between Is and 3s, — 
absorption is suppressed by the detuning A. 



Mathematically speaking, we assume that the Hilbert space of the atom is spanned by 
the three orthonormal basis states 

?/atom = hn{|ls),|2p),|3s)}. (35) 

When considering the well-established Hamiltonian for an electron bound by a potential, 

ff. = |1 + V(,) , (36) 

El, E2 and E^ are defined by the following characteristic equations 

i^olls) = ^i|ls) , i^o |2p) = ^2 |2p) , ifo |3s) = ^3 |3s) . (37) 

The level spacings are shortly denoted by £'12 = E2 — Ei, E23 = E^ — E2 and 
£^13 = E^ — El. In order to be able to observe two-photon absorption, we assume 
two incoming photons, e.g. the target photon with frequency Ui and the control photon 
with frequency 002- Both photons together are in resonance with the transition between 
Is and 3s, i.e. £'13 = wi + (X'2. Both photons by themselves however are detuned from the 
transition between Is and 2p, such that a one-photon absorption is forbidden by energy 
conservation. The same applies for one-photon transitions between Is and 3s, which are 
moreover forbidden by angular momentum conservation, as a photon always carries an 
angular momentum, but the involved states have the same angular momentum (/ = 0). 
To sum up, the level scheme supports two-photon absorption but inhibits one-photon 
absorption. However, one-photon scattering will remain as the dominant one-photon 
loss mechanism. 
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3.1. Hamiltonian and initial state 

We start most generally from the non-relativistic Hamiltonian for atom-field interactions 
[6T| {h = c = Eq = iiQ = 1), dividing it into the non-interacting part Hq and the 
interaction part Hi 

1 , 
2m 



H = ^(p-qA{r^ + 5r, t)) ' + V{r) + Afield = Ho + . {3t 
2m V / 



The non-interacting part Hq carries the energy of the electronic states as well as the 
field energy and is given by 

^0 = + El I Is) (ls| + E2 |2p) (2p| + ^3 |3s) (3s| . (39) 

k,\ 

The interaction Hamiltonian, 

2 

H^ = ^^p.A{ro + 6r, t) + |- A Vo + Sr, t) , (40) 
m 2m 

contains the rest mass m and the charge q = — |e| of an electron, together with its 

momentum (vector) operator p and the vector potential operator for the quantized 

electromagnetic field (here written in the interaction picture) 

A (ro + Sr, t) = J2 9k^k ^e-''"''+'''<''''+^^^ + H.c. , (41) 

fc,A 

at the position Vq + Sr. Here is the position of the atomic nucleus and Sr is the 
position of the electron relative to the nucleus position. The amplitude is given by 



9t = -hh7:rT^^4- (42) 

\j 2 (27r) ujk 

Hence the interaction Hamiltonian ( HOj) describes one-photon transitions between 
adjacent levels by its mixed term oc p ■ A (ro + Sr,t), and possible direct two-photon 

- 2 

absorption or a one-photon scattering by its quadratic term oc A (ro + Sr, t). 

The Hamiltonian f l40|) can be rewritten in terms of atomic transition operators 
(j^^ = |ls) (2p|, a^^ = |2p) (3s| and a^^ = |ls) (3s| and photonic annihilation/creation 
operators d, Jd\ , . 

fc,A' fe,A 

jyu(r,.«) = E (j?.*"e--l" + <,*^'e--?') al^e"-- + H.c. 

fe,A 
fc,A 

/ ^ ^kk,\\ k,\ k,\ 

k,k,X,>' 

+ V 0^- -e'(' 

k,k,X,\ 

+ i^[7(r'o,t). (43) 
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The derivation of the associated couphng constants gj^^ etc. can be found in 



Appendix B.l Note that the Hamiltonian fH3l) is written in the interaction picture 



and therefore now carries phases due to the atomic energies. The remaining part 



H[^[^{rQ,t) is given in Appendix B.l because the terms it contains are unimportant 



for the upcoming analysis. Abbrevations A^^ := E12 — oJk and A^^ := i?23 — are 
introduced to indicate the detunings between the photonic energies and the atomic 
transition energies. The interaction Hamiltonian fH3|) works on an incoming state which 



consists of the atom in the ground state |ls) and two incoming photons a\ and a 

IV-in) =aI4|0)|ls) . 

The two incident photons are generally smeared out in fc-space, 
XI /1/2 (^) «i 



it 



"1/2 



in such a way that they respect the usual commutation relations 



0^2, 0-2 



and 



0. 



In order to do so, /1/2 need to fulfill the conditions 



k 



\fy2{kf = l and 5^ / W (fc) / W (fc) = 



(44) 



(45) 



(46) 



(47) 



This can be seen by employing the fundamental commutation relation for creation and 
annihilation operators, 



6Uk-k]6 



(48) 



3.2. Perturbation theory 

The time evolution of the system for infinite interaction time is governed by 



If {vq, 00) = Texp ( — i 



dt'H^^i{ro,t'] 



(49) 



with time ordering operator T and atomic (nucleus) position Vq. We calculate 
the outgoing state ?7(ro, 00) \ipi^) via second-order standard perturbation theory [56] , 
dropping all terms of order A or higher, lipi) and |^/'2) represent the first-order term 
and the second-order term, respectively, such that 

f/(ro, 00) l^in) ^ l^in) + 1^1) + IV'2) . (50) 

For infinite interaction time, only terms matter which satisfy energy conservation, 
because other terms have a non- vanishing temporal phase which makes them cancel out 
in the infinite time integral. Considering this (using also Ui/2 < E13) and performing 



some commutations of the annihilation- and creation operators (see Appendix C.l ), we 
arrive at 



iV^i) = -i(27r) 2{gll J*6{u, + u;-,-E,,) 



X 



fe,fc,A,A 
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x/;(fc)/*(fc)5,^A.Ae'^'^'^"'|0) |3s)- 
-i(27r) 9llx,H^k - i^k) X 

X [fm^A + mk)s,,-A) ^I,/^'"'^-"" io) lis) . (51) 

That is, the first-order term lipi) contains one-photon scattering and two-photon 
absorption via the oc A {r^ + 6r,t)-teYm of the original Hamiltonian fHOl) . see also 
Figure [71 



UJ2 



Eo 



E3 
E2 



1^1/2 



El 



(a) 



El 



(b) 



Figure 7. First-order processes of the Hamiltonian (|43p acting on the input state 
„ 2 

(HJ). The absorption (a) or scattering (b) works via the cx A (tq + Sr,t)-term of the 
Hamihonian (|40l) . thus the the 2p-level is irrelevant here. A circle with black dashed 
border symbolizes the initial state of the atom, while a green circle symbolizes the final 
state. For scattering, (b), initial and final state of the atom are the same. 



The second-order term is also given by fH9l) . where we additionally have to attend for 

normal ordering of the photonic creation and annihilation operators [62]. As stated 

^ 3 

above, we leave out all terms of order A or higher, so for second-order perturbation 
theory, only the first two hues of the Hamiltonian (H3|) contribute. Calculating the time- 
integrals, using the commutation relations (Appendix C.l ) and renaming indices in the 
first term we arrive at 

\ij2) = 1 (27r) 2^ I + 1 H^k + (^k- ^13) X 

X /;(fc)/*(fc)5,^A.Ae'^'^'^-^" |0) |3s) + 

fc,fc,A,A \ « / 

X {m~k)S,,A + m~kK.x4) |0) lis) . (52) 

In the second-order term, we also have two-photon absorption and one-photon 
scattering, but this time going over the intermediate level 2p via the oc p ■ A -term of 
the Hamiltonian (110|) . Therefore the detuning of the two photons from the intermediate 
level, e.g. Aj^ = E12 — u^, also plays a role here. Note that there are two terms 
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for two-photon absorption (bracket in the first hne), as there are two possible orders 
of absorption. There are also two scattering mechanisms (bracket in the third line), 
whereof the second one has a detuning comparable to optical energies and thus will be 
neglected. All four processes are shown in Figure [HI 



0J2 




& 



IA12| 



(a) 



^3 




(b) 

Figure 8. Second-order processes of the Harailtonian (|43l) working on the input state 
(H^ . The absorption (a) or scattering (b), exemplarily shown for wi getting scattered) 
works via the oc p • A -term of the Hamiltonian pp]) . A white circle with black dashed 
border symbolizes the initial state of the atom, while a white circle with green solid 
border symbolizes the occupation of a virtual level (black dashed line) during the 
transition. A filled green circle denotes the final state. The detuning of the respective 
process is given by the distance between the 2p-level (black solid line) and the virtual 
level (black dashed line). 



3.3. Two-photon absorption probability 



Using both the first-order fl5T!) and the second-order (152!) term, we can calculate the 
total two-photon absorption probability P27 

2 



27 



(3s| f/(ro,oo) iV^in) 



(53) 



In order to obtain a concrete result we now need to specify the distribution of the 
incoming photons in fc-space. For the sake of simplicity, we assume rectangular functions 
around fci/2 = ki/2 ■ e^, i.e. the incoming photons propagate almost parallel to the z- 
axis, with slight deviations in x- and ?/-direction. The side lengths are denoted 2Akx, 
2/S.ky and 2Afc2, and of course /1/2 are normalized such that (H71) is fulfilled. 
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As the side lengths in fc-space should be small, the integrand in (1^ can be regarded 
as constant over the small integration volume at k ^ and k ^ k2. The level scheme 
of the atom is designed in such a way that the detuning of the target photon (A^^ =: A), 
is much smaller than the detuning of the control photon (A^^), whose term is therefore 
omitted. For small detuning, i.e. in the optical regime, the direct transition [g^^ ^ ^ )* , 

produced via the quadratic oc A -term, can also be neglected ( [Appendix B.lj ). The 
result is 



27 



(2vr) 



16 



,23 



A2 



-AktAkl 



(54) 



Inserting the expressions for the coupling constants (notably the typical coupling length 
•^atom = I (^outl ^ li^in) \ , scc Appendix B.l for details) and introducing the cross section 
area of the two photon wave-packets A := 1/ (Afc^AA;^) finally yields 



4"qED -£^12-^23 



p4. 

, _ ^--^gED -^12-^23 ^atom /p-p-N 
7T'^UiU2 A"' 

where, for simplicity, we assumed that the polarization vectors are along the same 
axis as the absorbing orbitals. Otherwise, the result is still a reasonable estimate for 
the magnitude of the two-photon absorption. Note that the probability fl55|) does not 
depend on the length of the photon wave-packets; only the transversal dimensions (i.e., 
A) enter. 



3.4- One-photon scattering probability 

From the final amphtudes (ISTIl and (J52l) , we also may calculate the one-photon scattering 
probability Pi^. 

2 



Pi 



17 



;is| f/(ro,oo) l^in) 



(56) 



Using the fundamental commutation relation, one can verify that Pi^ can be simplified 
to the following form due to the orthogonality in k- and A-space. The index ( = 1,2 
sums over the two possibly scattered photons. From the commutation relations there 
also arises a term of order which will be neglected, as it would be about a factor 
O {AkxAky/ (fci/c2)) smaller than the first term. As stated above, the second scattering 
mechanism was ignored due to its very large detuning E12 + uj^- 

,12 / 12 Y 



I7 



(2vr) 



fc,A,C=l,2 



9 



A12 

k 



9^ 



11 

fcfc,AA 



C 



i{k—k)-ro 



+ o if) 



(57) 



After carrying out the inner integration (by assuming the integrand to be constant at 
k ~ fc^), the outer integration can be done in spherical coordinates where the radial 
integration is restricted due to the Dirac delta function in the inner integration (for 
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further details, see Appendix C.2). After inserting the expressions for the couphng 
constants and executing the outer integration and the sum over the two polarization 
modes, we arrive at the final result 

2 



17 



QED 



E 

C=l,2 



-^12 



A 



12 



^^atom 



?4 
'atom 



A 



(5^ 



It can be seen, that the one-photon scattering probability is based on two mechanisms. 
The first term in the bracket corresponds to a scattering via a virtual occupation of the 
intermediate level 2p, while the second term corresponds to a direct scattering via the 
oc A -term of the Hamiltonian fHOj) . Usually, for optical energies Eu and very small 
detunings (but still much larger than the natural line width of the middle 2p- level), the 
first term is some orders of magnitude larger than the second term, which therefore can 
be neglected. The first term for ( = 2 describes the scattering probability of the control 
photon. Choosing the detuning of the control photon to be much larger than that of 
the target photon, i.e. 3> A^^ = A, the scattering of the control photon can be 
neglected as required in section [2. 2. 31 



3.5. Destructive interference 

In principle, it is possible to suppress the scattering by destructive interference between 
the two contributions (from the oc p ■ A- and the oc A -term, respectively) to the 
scattering amplitude, leading to Pi^ ^ P2-y. To this end, we need to go to the infrared 
regime and to large detunings, such that the first term in ( !58l) gets small. When dealing 
with large detunings, neglecting the second scattering term in ( l52l) is not accurate and 
therefore we keep it for this consideration. For the interference to be destructive we 
obtain 



12 / 12 



+ 



12 



9 



11 



0, 



(59) 



Whereas in the previous sections it was convenient to assume that the polarization 
vector of the incoming photon e^^ is along the same axis as the absorbing orbital 
{ls\r\2p), it is necessary for destructive interference to happen. So when 



pl2 
'atom 



0I2 



we require this and additionally l^^atomi 
E^ f ' ^ 

-^12^atom 



?atom, our condition fl59l) simplifies to 



1 



E- 



12 



UJkr 



E12 + UJkr 



pl2 

"atom 



''atom 



e 

m 



fee 



The expression inside the square brackets is zero when 



E- 



2 «2 
12^atom 



E 



12 



E12 + UJk. 



1 

= 0, 

m 



(60) 



(61) 



or written as a condition on the frequencies of the target photon and the control photon 
ui=U2 = Eu Vl - 2m£Lm^i2 . (62) 
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Together with the resonance condition UJ1+UJ2 = £"12 + £"23, there is only one parameter 
left. For example, choosing E12 determines Ui = UJ2 by fl62|) . Due to the resonance 
condition, £^23 is also determined then. 



3.6. Comparison and conclusion 

In the general case without destructive interference, one can neglect the second term 
inside the bracket of (158|) . as well as the first term for C = 2 due to the large detuning 
of the control photon. What remains is the scattering probability of the target photon 
due to the oc p ■ A -term of the Hamiltonian pUj) . which is 



^^^-^^A^^- ^^^^ 
We arrive at an expression for the quotient of the two-photon absorption probability 
and the one-photon scattering probability 



(64) 



To get an impression what this ratio is in view of the diffraction limit, we can insert 
w ?^ Wi ~ (X'2 ~ E12 ~ E23 (for simplicity) and A = (A/2)^ where A = 2tt/uj, which 
results in 



P27 3 



0.05. (65) 



So even at the diffraction limit, this ratio is smaller than one and thus the requirement 
f l28|) does not seem to be satisfiable without additional efforts. In the next sections HI 
m and O we present possible mechanisms to boost two-photon absorption in order to 
overcome this difficulty. 



3.7. Effective treatment via the "slow variables" formalism 

In our analysis of two-photon absorption, section l33| we found that the two-step process 
via the (x p ■ A -term is the dominant process. In what follows, we will thus repeatedly 
resort to an effective treatment of two-photon absorption and emission, where we omit 
the middle (2p) level from the Hamiltonian and model the two-step process effectively 
as direct transition between the Is- and the 3s- level, first-order process. 

The "slow variables" formalism, which we explain below, provides an elegant 
way to pass into the effective treatment and to evaluate the associated coupling 
constants. Considering our three-level system, we start from a reduced Hamiltonian 
which treats the electromagnetic field classically, supporting only the modes Ai and 
A2 with frequencies ui and U2, respectively. Moreover it supports only transitions 
between adjacent levels via the appropriate mode Ai or A2, i.e. the cx A -term is being 
completely disregarded 

Hred = El I Is) (ls| + E2 |2p) (2p| + E3 |3s) (3s| + 

+ {g^^Aia^\''''' + ^23^;a23gi^2* + H.c.) . (66) 
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From this reduced Hamiltonian, we can construct an associated Lagrangian which treats 
the state-vectors as if they were complex numbers (which are of course normahzed) 

+ {g'^Al^l^^pe''^'' + g''A;%^^sse''''' + H.c.) . (67) 

Redefining the phases of the complex variables in order to eliminate explicit phases in 
the interaction terms and to cancel out two of the three energy terms 

V'ls = ^ise^""^*, V'2p = *2pe'(^^+'^^)*, ^3s = *3se'^^* , (68) 

yields a shorter Lagrangian 

L = A^/'^p^ap - i (^LV'ls + ''/'2p^2p + ^3s^3s) + 

+ {g'^Airis^2p + 9''A;r2p^ss + H.c.) , (69) 
which involves the detuning A = E12 —Ui. From (IMj) the Euler-Lagrange equations for 
the variables ip2py i^is V'sp '^ead 

A^2p - i^2p + g'^*A,ij,, + g^^AliJss = 0, (70) 
-iij,, + g'''Al^2p =0, (71) 

-i^/^ss + g^^*A2iJ2p = . (72) 

fITO]) can be solved approximately for ilj2p in the following way 
^P2p = - (A - idt)-' {g'''*A^^u + 9'''A;^Pss) 
= - ^ {9''*A,4j,, + g''A;^p,,) - 

- ^ {\9''*A,\' + \g''A;\') 4^2, + O ilgAf/A') , (73) 

where the Euler-Lagrange equations for ipl^ and ip^^ have been inserted. We assume 
that \gA\ <^ A and thus 

V^2p ^ - ^ (^?''*Ai^is + 9''A;^,,) . (74) 

Reinserting fTMl) into the Lagrangian fl69l) . where the '?/'2p-term was neglected as well for 
the same reason, yields an effective Lagrangian in which the 2p-level was integrated out 



LcS= - i (V^LV'ls + ^^38^^38) 



\g'^*A,ij,, + g^^A*iJssr 
A 



/ 12* 23* \ 
+ f AiA2^3*,V^is + H.C.J . (75) 

We find the two-photon absorption and emission processes in the last line. The coupling 
constant of these processes is thus given by 
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4. Repeated inducing 

In section [21 we presented our two-photon entangling sclieme whicli requires the two- 
photon absorption rate of the absorbers to be much larger than the one-photon loss rate. 
We subsequently investigated the situation in case of a simple three-level atom (section 
[3]), and obtained the result that the two-photon absorption probability is at best still 
about 20-times smaller than the one-photon scattering probability. Therefore, additional 
enhancement of two-photon absorption is necessary for our scheme to work. Thus in 
this section, a mechanism is presented which could enhance the two-photon absorption 
probability compared to the one-photon scattering probability. The mechanism is again 




Figure 9. Internal structure of the absorbers where two photons are sent through the 
same absorbing medium repeatedly (in this case, n = 4 times). The photons always 
travel the same length L between two successive passages of the absorbing medium, 
which was disregarded in this picture. In order to ensure that the incident angle of the 
photons remains the same, the medium could be rotated. 



based on the quantum Zeno effect. When sending the two photons n times through the 
same absorbing medium (see Figure [9]) with the correct optical path length L, the two- 
photon absorption amplitudes add up coherently, while for scattering processes, only 
the probabilities add up. To see this, let us first introduce an effective Hamiltonian 

V;!?^ = 9i3CT^^alal + gi2a^^a\ + fea^'a^ + giMcluer) (V-inl + H.c. . (77) 

In this effective Hamiltonian, for the sake of completeness, we also included one-photon 
absorption and emission (the case of zero detuning). The index denotes that this 
is the Hamiltonian for the first passage of the absorbing medium. The incoming state 
is ( l44l) as before. In the outgoing state, there are small amplitudes for two-photon 
absorption, one-photon absorption and one-photon scattering according to ( 177]) . We 
assume small interaction time r <C 1, i.e. we regard only up to first order in r 



= iV^in) - ir [gl, |0) \3s) + gl.&l |0) \2p) + ^?n|V^L°Lr)J + O (r)^ . (78) 

Afterwards, the two photons are guided into the material again, crossing a certain 
optical path length L, which should be much larger than the size of the photon wave- 
packets. The Hamiltonian for the second passage, v}^\ therefore carries different phase 
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factors depending on which photons are involved in the respective processes. That is 
because the Hamihonian for the second passage couples to the vector potential at some 
position A(ro + rL,t) instead of A{ro,t), and the spatial phase factors expjifc ■ vq} 
change accordingly. In contrast to the two-photon absorption, the scattering process 
behaves incoherently, i.e., like repeated measurements, when L is larger than the size 
of the photon wave-packets. There only the probabilities add up and the scattering 
amplitudes of the different passages are orthogonal. Hence they are labeled with indices 

(V'scatterlV'scatter) — ^ij 

+ 9l2e-"''''a'^ + (?23e-''^^a'^a^ + 9u\4cLor) {tj + H.c. . (79) 

So in the second passage, the state gathers additional amplitudes with their respective 
phase-factors 



irV;W + O (r 



= iVin) - ir gl, (1 + ei('=i+'=^)^) |0) \3s) + gl, (l + e'^^^^) a| |0) \2p) + 



+fll| Ater) + ^lllCatter) 



(1) 



+ 0(r) 



(80) 



When the two photons are repeatedly guided into the absorbing material for n times, 
each time passing the same optical path length L, the state evolves as 

n-1 



n-1 



= iV'in) -irJ2 [^i*3e^('^+'^)'^^ |0) \3s) + 

fM=0 

+ gl^e'^'^'^dllO) |2p)+5n|V^ilter 



+ 0(r) 



(81) 



Moreover, when calculating the overall two-photon absorption probability, we readily 
see that it is enhanced by a factor n^, assumed we choose the optical path length such 
that is related to the sum of the two photon wave-numbers via {ki + k2)L e 27rN 



27 



|(3s|^(")) 



2 I * I 



n-1 



2^2 



(82) 



A possible local one-photon absorption effect, in contrast, would violate the phase 
matching requirements when choosing kiL ^ 27rN. Therefore the one-photon absorption 
probability does not even scale with n. It is heavily suppressed 

2 



pabs 
^7 



|(2p|*W)f ^T'\gl,\ 



n—l 



^ \9l2\ 



^l=o 

2 cos nkiL — 1 
cos kiL — 1 



1^1*2 I' 0(1). 



(83) 
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Most important, for one-photon scattering only the probabihties add up 



n-l 2 



pscatter ^ ^ | (a')^^^ J vl/ W ) = \gn\' n . (84) 

At=0 

In this way, we can enhance the total two-photon absorption probability by a factor of 
as compared to the one-photon scattering losses, which scale with n. Concluding this 
section, we see that the requirement (!28|) could be achieved by sufficiently large n 



67 _ n 
5. Coherent excitation of many atoms 

It was already pointed out that additional efforts are needed to further enhance 
two-photon absorption as compared to one-photon scattering. One mechanism to 
achieve such an enhancement was already given in section HI In this section, we 
present an alternative approach based on coherent excitation of a large number 5* of 
atoms/molecules. We will show that the two-photon absorption probability scales with 
the number of excitations, whereas the one-photon processes do not. In addition we 
discuss how an excited state of many atoms/molecules could be sustained by pump 
lasers. Coherently excited states of many atoms are often referred to as "Dicke- 
states" [63], which are typically investigated regarding collective spontaneous emission 
also known as "Dicke super-radiance" [64] - [72] . As we now deal with many three- level 
atoms instead of one, our Hamiltonian is now the sum of the individual single-atom 
Hamiltonians in (HSl), i.e. 



^ciustcr({ra, t) = J2 HiAfi. t) . (86) 

e=i 

The single atomic transition operators inside the individual single-atom Hamiltonians 
Hi^i{r£,t) are now denoted with an index i, e.g. aj^, which means that they act on the 
i-th atom. The initial state changes in so far, that we start from an 5'-atom state \s) 
instead of the ground state of a single atom |ls), 

\^l;t,) = a\al\0)\s) . (87) 

The s in |s) means that out of S atoms, s atoms should be coherently excited to the 
3s-level in order to enhance two-photon absorption, while S — s atoms are still in the 
ground state |ls). None of the atoms in \s) is in the state |2p), which corresponds to 
the middle (2p) level. 

5.1. Effective Hamiltonian and quasispin operators 

To get the idea behind coherent excitation, it is sufficient to look at the following 
effective (interaction) Hamiltonian, considering S atoms/molecules with their respective 
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positions 

-^Dicke = 5'«i«2 ^ ^l^'' exp{ir^ ■ (fci + ^2)} + H.c. . (88) 
It can be understood as reduced version of ( 186|) . with the middle (2p) level integrated 



out (resulting in an effective coupling constant g, see Appendix B.2 ) and all effects 
disregarded except (coherent) two-photon absorption and emission. As we only have 
the Is- and the 3s-level left in this effective treatment, the Hamiltonian ( l88|) describes 5* 
two-level systems. These S two-level systems can also be regarded as S spin-| systems, 
identifying the |ls)-state with spin down, = — |, and the |3s) state with spin up, 
Sz = Accordingly, the atomic transition operators correspond to ladder operators 
from Pauli matrices af = {a^ ± i(3"|)/2, such that a/ = a^^"^ and = a}^. With this 
in mind, we can define quasispin-S' operators which sum over the S individual spin-| 
systems as 

s 

±x ^ .j.y ^f^± af exp{±ir^ • (fci + fca)} . (89) 

£=1 

Moreover, for the S spin-i systems, we also have the af-operator from the Pauli-z- 
matrix, with eigenvalues ±1, such that (t| |ls)^ = — |ls)^ and (5"| |3s)^ = -|- |3s)^. We 
define the quasispin-S* operator 
1 ^ 

S-=-^a|, (90) 

1=1 

which essentially sums up the energy due to excitations, Hq = i?i3(S^ + S/2). The 
Hamiltonian fl88|) can then be written shortly as 

-^Dicke = gaid2t^ + H.c. . (91) 

The operators S^, ±y and ±' form an SU{2) algebra \±^,±''] = Epi^^^^P^^ PS]- 
Therefore they have the same characteristics as the usual spin operators and especially 
they have the same normalized eigenstates. As is well-known [56], the following relations 
for the ladder operators arise solely from the angular momentum operator algebra 

X]p ^^fiupjp- The |j, m) form a set of normalized eigenstates of J and 
J+ 1 J, = - m){j + m + l)\j, m + 1) , 

J_|j, m) = a/ (j + m){j -m + m - 1) , (92) 

where the quantum number m can take values from —j to j with unit steps. 

The same is true for our quasispin operators S^, and S^. The variables however 
are defined slightly different. That is, the total number of the individual spin-^ systems 
S corresponds to the maximum value of m, i.e. S = 2j, and the number of excitons s is 
related to the quantum number mhjs = m + j,m ranging in —S/2 to +5/2, s takes 
values from to S. This results in the following normalization factors 

±+\s) = y/(S -s){s + l)\s + l) , (93) 

t-\s) = y^iS -s + l)s\s-l) , (94) 
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where the state |s + 1) denotes an entangled state of s + 1 coherently excited atoms and 
S — s — 1 atoms in the ground state. Analogously, |s — 1) represents a state where s — 1 
atoms are coherently excited to the 3s-level and — s + 1 atoms are in the ground state. 



5.2. Enhanced two-photon absorption probability 

Regarding particularly ( !93ll . one sees that the absorption probability scales with Ss in 
the limit s <^ S* 

P^^ = {S-s){s + l)P2^^SsP2^. (95) 

This means that the two-photon absorption probability is enhanced by a factor s, as it 
would only scale with S in the usual case without coherent excitation, i.e. s = 0. Please 
note that the corresponding directed two-photon emission (see [65l[67l[7l]) is amplified 
too, (Elj), but this means no harm to our desired two-photon absorption, as long as the 
atoms can be kept in the excited state, e.g. by constantly pumping (see below). 

Aside from two-photon absorption and emission and one-photon emission processes 
(see above), we find one-photon scattering on the ground state as well as on the excited 
state, and a special upconversion process where one photon is absorbed, one atom relaxes 
into the ground state, and a higher energetic photon is emitted again. As for the single- 
atom case, full perturbation theory calculations have been performed for fl86|) in order 
to maintain the probability for any possible process. All important processes are shown 
in Figure [TOl 



5.3. One-photon scattering probability 



Let us now take a closer look at the one-photon scattering processes. The amplitude for 
the one-photon scattering on the ground state is proportional (taking only the atomic 
part, the rest is comparable to the scattering amplitude in section [3]) to 

s 



(96) 



where k is the fc- vector of the scattered photon before the scattering, e.g. k = ki, and 
k is the wave vector after the scattering. For k ki, the phases inside the sum are 
randomly distributed because of the random atomic positions r^. Therefore the phases 
do not add up coherently but rather statistically (like in a two-dimensional random 
walk) 

2 



^exp{ir£ • (fci - k)} 



eel 



0(3) . 



(97) 



/ denotes an arbitrary index set of {S — s) atoms in the ground state. Thus for the 
scattering probability in case of S atoms we simply arrive at 



l7 



O (S) P, 



I7 5 



(9^ 
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Figure 10. All important processes (coming from the p ■ A -term in the Hamiltonian) 
in the situation of a cluster of coherently excited atoms |s) are shown, (a) Two- 
photon absorption, probability cx Ss. (b) One-photon scattering on the ground 
state (exemplarily shown for uji getting scattered), probability oc 0(5*). (c) One- 
photon scattering on the excited state (for 0-12), probability cx 0(s). (d) One-photon 
upconversion (for wi), probability oc Ss, but very large detuning. 
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where Pi^ was the scattering probabihty for one atom. The same argument apphes for 
scattering on the excited state, except that there are only s atoms in / and thus the 
scattering probability is only oc O (s) Pi-y. 

If, however, k = ki, (!97|) evaluates to {S — s)^ and the scattering seems to be 
heavily enhanced. But in this case the photon is exactly in the same state as before, i.e. 
it actually has not been scattered but just passed through the atomic cluster. In other 
words, nothing happened. 

To summarize, the one-photon scattering processes are at most enhanced by a factor 
0(5") compared with the single-atom case (section [3]). This result is also intuitively 
plausible, as the photon can now be scattered by 5* different atoms instead of just one. 
Notably, the scattering probability does not scale with the number of excitations s, as 
the two-photon absorption, since it does not satisfy the phase matching condition. 

5.4- One-photon upconversion probability 

In this process, a single photon gets collectively (virtually) absorbed by the excited 
atoms, followed by the emission of a higher-energetic photon via relaxation into the 
next lower coherent state |s — 1). The emission is also directed due to the collective 
character of the transition, i.e. the phase matching condition. For example the target 
photon ki is (virtually) absorbed, followed by the emission of a photon with wave 
vector 2ki + k2. This can be seen by looking at the atomic part of the amplitude of this 
upconversion process 

s 

l^conv) oc J2 ^M-^re -{k-k)} \s) . (99) 

Again, k is the fc- vector of the photon before the absorption. Sticking to our example, 
k = ki, the scattering is directed along k = 2ki-\-k2 in order to fulfill the temporal phase 
matching condition k — k = ki + In that case, the atomic part of iV'conv) given by 
\iplonv) OC S~ |s) = ^y {S — s + l)s |s — 1), and therefore the probability for this process 
also scales with Ss in the limit s S. Fortunately, this process is not supported by the 
level scheme of the atom. In other words, the detuning for this process is comparable to 
optical energies. From Figure fTOT d) it can be seen that for our example, the detuning is 
E23 + Ui respectively —E12 — Ui (in the case where the control photon U2 gets absorbed, 
it would just be E23 + uj2 and —E12 — uj2)- Therefore, this process does not scale with 
1 / as in f l55|) and thus can be neglected for small detunings. 

5.5. Comparison and conclusion for many atoms 

Summarizing, the two-photon absorption probability is enhanced by a factor s compared 
to the one-photon scattering probability 
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Concluding this section, we see that the requirement (^E^ could be achieved by 
sufficiently large s 

K = i^ = -— -i^. (101) 
5.6. Sustaining the excited state 

As discussed above, the state of coherent excitation of s atoms, decays very rapidly 
by directed, spontaneous two-photon emission. In order to sustain the state |s), we 
propose to apply two pump lasers with wave numbers fci', k-^ which satisfy the same 
spatial and temporal phase matching conditions k\' + = ki + k2 and ijj[+U2 = ijJi+co2 
as the control and target photon. The quotient of the atoms which can be kept excited, 
s/S, should be primarily depending on the strength of the pump laser fields, which are 
naturally limited. To obtain a relation between the quotient s/S and the intensities 
of the pump lasers, Ji and I2, we can start from a Hamiltonian similar to ( lOTl) . but 
regarding the pump laser fields classically. Moreover we write the full Hamiltonian, 
including the energy due to excitations, and not only the interaction part 

i^Laser = ^13 + + (^gt+ + H.C.) . (102) 

In this case, the coupling constant g reads (see [76]) 

1,1 . ^J>^A,A, = *Jl^^SE^^^lA^A,A, , (103) 

where A' = E12 — u)[ denotes the detuning of the pump laser field. We can then treat 
the quasispin-S* system as a harmonic oscillator by applying a Holstein-Primakoff [71] 
transformation, 

t 



S+ = a^^S-a^a = (t-J , t' = a^a-S/2, (104) 

on ffT02D . This yields 

^Laser = ^i3a"^a + (^gaW S - ata + H.c.j . (105) 

In our envisaged limit S ^ s the Hamiltonian fllOSp looks just like a harmonic oscillator 
coupled to the classical laser fields 

^Laser = ^13«^« + {ga^ + H.C.) . (106) 

Given a certain coupling constant g, which particularly depends on the laser- field 
intensities and the detuning, there is an associated coherent, steady state d \ag) = 
ctg |ag), which can be found by solving the characteristic equation HLaser l^g) = Ea |ag) 

-^Laser 1%) = (^-ElsCtg + V^5') «^ + '/Sg*^^ jOg) = E^ l^g) . (107) 

Thus we found the eigenstate of the S'-atom system, it is the coherent state with 

ag = -^v^/Ei3, (108) 

and 

E„ = -Ei3|«g|'. (109) 
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The average exciton number in a coherent state is given by (n) = {a'^d) = |agp. Thus, 
by setting |ag| = y/s, we get an expression for the quotient s/S 

The vector potential is related to the intensity by / = A^u"^ (remember c = eo = 1) and 
we may express E13 by the two frequencies E13 = uj[ + 002- Thus the intensities of the 
two pump lasers, Ji, I2, are related to the ratio s/S via 

S - [ co'M + u;^) A' ) ■ ^^^^^ 
The detuning A' of the pump laser field cannot be chosen too small, because this would 
lead to non-negligible population of the middle (2p) level which we integrated out in 
our Hamiltonian ( I102p . The population of the intermediate level is governed by 

V'2p = {g''*A,ij,, + g'^A^iJss) . (112) 

In other words, the detuning should be big enough such that g^'^*Ai/A' ^ 1 as well as 
^23^*^^/ ^ j g nearly no population of the middle (2p) level. In our Hamiltonian 
notation this corresponds to \g^'^ \ Ai ^ A' and |gr^^| A2 ^ A'. Inserting \g^'^\, \g'^^\ and 
again / = A'^u'^ gives 



V'47r«QED/l |4'om| E12/UJ1 < A' , 

V47raQED/2 l^oml ^23M « A' . (113) 

Assuming that Ei2/u)i ~ i?23/i^2 ~ 1 and that the laser intensities for both pump lasers 
are the same, Ii = I2 = I, we find that the detuning A' of the pump laser field should 
satisfy the condition 

atom ; 

(114) 

in order to avoid unwanted excitations of the middle (2p) level. With a typical dipole 
length of six Bohr radii ^atom = Sob and for a large but possibly realistic intensity of 
/ = IQiow/cm^, this translates into A' > 10^^ Hz or A' > 0.06 eV. 



6. Example values 

Let us insert some typical parameters. We assume that target and control photon are in 
the optical regime (say around 500 nm) and adjust the detuning of the target photon to 
be A = 3 ■ 10^^ Hz, which is several orders of magnitude larger than the typical natural 
line width of the middle (2p) level (in the range of GHz). Choosing the detuning of 
the control photon to be an order of magnitude larger, i.e., 3 ■ 10^^ Hz, the loss rate 
of the control photon can be neglected. We consider three different values of the error 
probability Perror, namely Perror = 0.5, Perror = 0.25 and Perror = 0.1. At first, we 
determine the required ratio of the two-photon absorption and the one-photon loss rate 
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K = ^27/^17 and the needed number of segments of the gate. Recall that in section fI72\ 
we already derived the relationship 



K 



^27 
^17 



2PL 



> 1 



(115) 



which, however, is only valid for ^ 1 and ^27 ^ ^17- For experimental realization, 
small values for and k are desirable. Thus we calculate some possible values for A^ 
and K numerically. As there is some margin for trade-off between A^ and k, we present 
three different possible choices of A^ and k, a choice with small A^, a balanced one and 
a choice with small k. For each possible realization of the set-up in Figure [5l we also 



p 

error 


A^ 




p 

error 


A^ 


n 


p 

^ error 


N 




50% 


8 


22 


25% 


20 


120 


10% 


50 


1 430 




10 


12 




25 


76 




60 


760 




40 


8 




70 


55 




160 


440 



Table 1. Example values of the number N of segments and the corresponding ratio 
K. For each value of the error probability Porror, three possible choices of N and k are 
given. For instance, to reach an error threshold of Perror = 50%, it would be possible 
to build a gate with = 10 and k = 12. However, when k = 12 poses too hard a 
challenge, it is e.g. also possible to go with = 40 and k = 8. 



specify a two-photon absorption probability P27^™ and a one-photon loss probability 
nt, which is related to the rat 

(see section 12. 2p 



psegm segment, which is related to the ratio k and the number A^ of segments by 



psegm 

^27 



psegm 
^17 



1 



exp 



exp 




;ii6) 



Furthermore, we state what is the number n of repetitions necessary to reach the 
respective ratio k via the mechanism described in section H] (sketched in Figure |9]). 
The required amplification n is simply given by (see above, section 13. 6p 

n^ = K. (117) 

-T I7 

This expression was evaluated using fl55|) and fl58l) . inserting the values given in this 
section and assuming a focus at the diffraction limit, A = (A/2)^. As the required 
amplification of the two-photon absorption probability can also be achieved via coherent 
excitation of a large number s of atoms/molecules (see section E]), the obtained value in 
(11171) can also correspond to s or even to ns, when both enhancement mechanisms are 
combined. The discussed values are given below in Table [2] for the choice of small A^, in 
Table [3] for the balanced choice and in Table H] for the choice of small k. To see which 
values of s might be realistic, let us insert a detuning of the pump beam of A' = 3-10^^ Hz 
into expression flllll) . where we get an excitation ratio of s/S = 1.7 ■ 10^^. In order to 
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^ error 


N 


psegm 


psegm 
^7 


t — 67/ 67 


n, s or ns 


50% 


8 


99% 


21% 


22 


471 


25% 


20 


99% 


4% 


120 


2 567 


10% 


50 


99.9% 


0.5% 


1 430 


30 588 



Table 2. Example values for the choice of small N. Listed is the number of segments 
N, the two-photon absorption probability P27^™ and the one-photon loss probability 
P^^^™ per segment, the corresponding ratio k, as well as the enhancement factor n, s 
or ns necessary for reaching this ratio. 



p 

error 


N 


psegm 
^27 


psegm 


— 67/ 67 


n, s or ns 


50% 


10 


95% 


23% 


12 


257 


25% 


25 


95% 


4% 


76 


1 626 


10% 


60 


98% 


0.5% 


760 


16 256 



Table 3. Example values for the balanced choice. Listed is the number of segments 
N, the two-photon absorption probability Pa^^™ and the one-photon loss probability 
^1°^'° per segment, the corresponding ratio k, as well as the enhancement factor n, s 
or ns necessary for reaching this ratio. 



p 

error 


N 


psegm 
^27 


pscgm 


— 67/67 


n, s or ns 


50% 


40 


47% 


8% 


8 


171 


25% 


70 


61% 


2% 


55 


1 176 


10% 


160 


69% 


0.3% 


440 


9 412 



Table 4. Example values for the choice of small k. Listed is the number of segments 
N, the two-photon absorption probability P^^^'" and the one-photon loss probability 
per segment, the corresponding ratio n, as well as the enhancement factor n, s 
or ns necessary for reaching this ratio. 



obtain a reasonable value for S, we imagine a glass plate of = 10 /im thickness, for 
example, and a cross section area of A = (A/2)^, where we assume A = 500 nm. With 
a density of Si02 of p ~ 2.5 g/cm^ and a molar mass of M = 60.1 g/mol [75], we find 

iVmolecules ~ ~ 1.6 ■ 10^° , (118) 

where Na is the Avogadro constant. If one percent of these molecules is optically active, 
we get 5* = 1.6 ■ 10*^ and thus s = 2 720. Inserting the aforementioned values for the 
detuning of the target photon (A = 3 ■ 10^^ Hz) and the control photon (3 ■ 10^^ Hz), 
the two-photon absorption probability in ( l55l) would be around P27 = O(10~^^), which 
roughly fits to the value of 5" = 1.6 ■ 10^ discussed above. 

Going to the limit, we may imagine increasing these two detunings to A = 3-10^^ Hz 
for the target photon and for the control photon 3 ■ 10^^ Hz (which is in the infra-red 
region). In this case, the two-photon absorption probability in ( l55l) is two orders of 
magnitude lower and we could use a glass plate of 1 mm thickness, which increases the 
maximum number s of excited atoms/molecules by two orders of magnitude. The values 
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in the three Tables remain the same with the exception of the last column, where the 
required numbers for n, s or ns increase by 25%. 

For the amplification mechanism sketched in Figure [9l a very tight focus is desirable. 
For the other enhancement mechanism based on (IHHIl . however, the spatial phase 
matching conditions become problematic if we focus down to the diffraction limit 
A = (A/2)^ due to the uncertainty in k. Therefore, let us consider increasing the 
cross section area A. On the one hand, this would decrease the ratio flMj) even further - 
but, on the other hand, the number S of atoms/molecules within this area A grows by 
the same factor. If we keep the pump laser intensity constant, the enhancement factor 
s compensates the shrinking ratio (Ell), i.e., a weaker focus with smaller uncertainty in 
k is also feasible. 



7. Alternative level scheme 



As another idea for suppressing one-photon loss in comparison to two-photon absorption, 
let us replace the level scheme discussed in sect ion |3] (see Figure [6]) by a three- level system 
as in Figure [TT] where the middle level \tp2) lives much longer than the upper level {ips)- 
Such a level scheme is often referred to as A-system. 




Figure 11. Sketch (not to scale) of the level scheme. Both photons together are in 
resonance with the transition between \^pl) and |'03): E^ — Ei = lui+u;2, but one-photon 
absorption is suppressed by the detuning A, where E2 — Ei = oji — A. 



In this case, the coupling strength gi2 of the {ipi) — )■ \ip2) transition is much smaller 
than the other coupling strength (723 of the \ip2) I'^s) transition. According to Eq. ( |5^ . 
the two-photon absorption probability P27 scales with gl2- In contrast, as shown in 
Eq. flSTjl . the one-photon scattering probability Pi^ behaves as gf2- Thus, if we denote 
the ratio of the two coupling strengths by / = gi2/g23 ^ 1, we see that k = P27/-P17 
scales with l/P, and consequently can be enhanced strongly. 

Let us insert some example values. We assume that the life-time of the upper 
\1Jj3) level is of the order 10~^s, which is a usual value for fast optical transitions. The 
life-time of the middle \ip2) level is supposed to be much longer, about 10~^s. This 
corresponds to a ratio of / = (712/(723 = 0.03. Note that such differences are not unusual 
in quantum optics. For example, the couphng constants of an electric dipole (El) 
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transition (with Al = ±1, cf. [76j) and an electric quadrupole (E2) transition (with 
Al = 0, ±2) differ by a factor of about / = ^atom/A (cf. [77]) which is even smaller 
than the value / = gi2/g23 = 0.03 we have chosen. Inserting this value, we get at the 
diffraction limit 

i.e., an enhancement by a factor of 10^ in comparison to ( l65|) . This value of k corresponds 
to an error probability as low as Perror ~ 25 % (see Tabled]). 

Of course one needs to keep in mind that the two-photon absorption probability 
(133]) is also reduced by a factor of = 9 ■ 10^'^ in comparison with the example values 
in section [61 This poses constraints on the number of A-systems required for achieving 
a large enough two-photon absorption probability. In principle, this probability can be 
increased by reducing the detuning A, but this reduction is limited by the line- width 
of the middle \ip2) level. Unfortunately, this is not the natural line-width given by the 
inverse of the life-time, but gets broadened, e.g., by inhomogeneities and thermal effects. 
Here, we assume a detuning of A = 3 ■ 10^ Hz (for the target-photon). This value is a bit 
lower than in that of the gain media in typical solid-state (e.g., ruby) lasers, but could 
be achieved with gas lasers (see, e.g., [78]) or ions in crystals [791481] . Inserting this 
detuning, the two-photon absorption probability is about P27 = 0(10^^) per A-system. 
An absorption probability of order unity can then be achieved by packing 10'' or more 
of these A-systems in the focal point with a volume of order 1 /im^. This corresponds 
to a density of lO^^cm"^ which is quite feasible. 

Since it is probably hard to increase the spatial density of A-systems much more, 
higher values of / would require smaller detunings A < 3 ■ 10^ Hz in order to maintain 
P27 = 0(1). As one possibility, one could imagine using ultra-cold atoms or molecules 
as absorbed media. Although their spatial density is also limited, the line-widths are 
typically much sharper. For example, lowering the detuning one additional order of 
magnitude (A = 3 ■ 10^ Hz) would allow for / = 3 ■ 10~^ and thus for k, = 5.4 ■ 10^ which 
would correspond to an error probability far below 10%. 



8. Summary and conclusion 

In summary, we presented a scheme for entangling photons via the quantum Zeno 
effect which might be realizable with present-day technology. Compared to the 
previously proposed Franson gate [50H5lj . our three-branch gate (section [221) exhibits 
a significantly reduced error probability (section 12.2.21) due to its improved design. Our 
scheme could be physically implemented using strong two-photon absorption in optical 
fibres, as already envisaged for the Franson-gate [57H60] . In this case, the benefit of our 
scheme is given simply by the reduced error probability. 

Moreover, we explored three mechanisms to enhance the two-photon absorption 
rate compared to the one-photon loss rate. Because wave-guides, resonators or optical 
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fibres also have their drawbacks (as they, for example, can induce additional losses or 
decoherence) , we mainly focus on a free-space realization of our entangling gate. 

First (section H]), we presented an apparatus (Figure [9]) which allows to enhance the 
two-photon absorption probability compared to the one-photon scattering probability 
by a factor n by sending the two photons n times through the same, given absorber. 
The whole apparatus then replaces the absorbers in our three-branch gate (brown circles 
in Figure [5]). Second (section [5]), we proposed to employ coherent excitation of a large 
number of atoms/molecules, as known as "Dicke super-radiance" [65lin3[ZT]. It was 
shown that the two-photon absorption probability is enhanced by a factor Ss where S 
is the number of atoms and s is the number of excitations. Furthermore it was shown 
that any important one-photon process only scales with the number of atoms S. Thus we 
achieve the desired effect, the two-photon absorption is enhanced by a factor s compared 
to one-photon processes. Third, we consider a special level scheme (A-system) in section 
[TJ which possesses a strongly reduced one-photon scattering rate. 

Finally, we inserted potentially realistic example values in sections |6] and [3, showing 
that the presented scheme is not out of reach experimentally. Even though the example 
parameters indicate that it might be hard to directly reach the error threshold of 10~^ or 
10~^ required for universal quantum computation [8], the achievable success probabilities 
for entangling photons are already comparably high. To prevent any misunderstanding, 
it should at this point be stressed that in this work, the (pseudo) deterministic creation 
of entanglement between two given photons was considered, which is very different from 
the spontaneous creation of entangled photon pairs via parametric down conversion. 

As an outlook for future developments it should be stated that the design of our 
set-up in Figure [5] does not raise the claim to be optimal. It is possible that improved 
set-ups can be developed which yield an even lower error probability, given the same 
ratio of two-photon absorption compared to one-photon loss. Another key point where 
improvements could ensue is by providing new ways to enhance two-photon absorption 
compared to one-photon loss. With regard to the mechanisms proposed in this paper, 
better mirrors (i.e. more repetitions possible, see section Hj), stronger lasers [to improve 
the ratio s/S, see fill II) ]. or absorber media with sharper line-widths (in order to reduce 
the detuning, cf. section [7]) are desirable. 

In order to explore another direction in the multi-dimensional parameter space, one 
could consider realizing the three-level systems via Rydberg atoms or artificial atoms 
such as quantum dots, which generate electronic bound states with a level structure 
similar to Figure [61 In this case, the dipole length £atom could be increased by several 
orders of magnitude, see ( lllip . but the realistic energy and intensity scales would 
also have to be modified. As another point, by placing photon detectors around the 
apparatus, one can detect occurring errors with a certain probability and thus apply 
error correcting schemes [53llS]. For example, Myers and Gilchrist demonstrated that 
the the performance of a quantum Zeno gate can be greatly enhanced using photon loss 
codes l53l . 
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Appendix A. Success probabilities of the optical set-ups 

Appendix A.l. Success probabilities of the two-branch gate 

The error probabilities for both operating modes of the two-branch gate, as defined in 
section 12.1.11 

2 



pl7 
erro 

p27 



(0, 1) ■ ^ol 



(A.l) 



are analytically given by ( IT2l) . and were expanded for the limit A ^ 1 in zeroth order, 
see (fT7|) and (fT8|) . When expanding them until first non- vanishing order in 1/A while 
keeping the assumption A,^27 ^ 1 ^ ^^^17, they read 

PeLr = N^lj + 67 + , (A.2) 



and 
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+ 



24A 



192A3e2^ 
+ O , 



+ 



(A.3) 



where in the two-photon case the first order in 1/A vanishes. 



Appendix A.2. Success probabilities of the three-branch gate 



The error probabilities for both operating modes of the three-branch gate, as defined in 
section 12.2.11 

2 



pl7 

erro 

p27 



(o,o,i)-v^, 



out 



:i,o,o)-v^, 



out 



(A.4) 



are analytically given by ( 123|) . and were expanded for the limit A ^ 1 in zeroth order, 
see fl26|) . Using Mathematica, they can also be expanded until first non- vanishing order 
in 1/A while keeping the assumption A^27 ^ 1 ^ ^^^17- Then they read 



i^JrLr = ^ + 67 + O (A-^) 



(A.5) 
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and 



p27 



TT 



+ 



47r2 - 371^ 27r^ + tt^ 
+ 



N^2y ■ Asm 24Ar3^27 



+ 



UN 

where again in the two-photon case the first order in vanishes. 



(A.6) 



Appendix B. Details on the coupling constants 

Appendix B.l. Derivation of the coupling constants 

In order to examine the corresponding couphng constants, we need to calculate the 
atomic transition matrix elements [77J. For transitions which are given by the mixed 
term oc p ■ A^Vq + 6r,t), dipole approximation is applied, i.e. p ■ A (rg + Sr,t) ^ 
p A [vq, t), and thus only the p-operator works on the electron states. For the transition 
between Is and 2p, for example, the transition matrix element g^^ reads 



g^^ = - — (Isl p |2p) = -— (Isl im 
m m 

igi?12^atom 5 



|2p) 



= iqEi2 (ls| r |2p) 
where the commutator relation [p^, f] 



-2ip was used and £ 



12 

atom 



(B.l) 
(ls| r |2p) is a 



vector with the typical dipole (coupling) length of about six Bohr radii iHom = 6ao and 
a direction depending on the orientation of the electronic orbitals Is and 2p. Taking 
also the amplitude g^* from the vector potential into account we arrive at a coupling 
constant 



\ * 



9 



12 



E 



12 



and analogously for g 



23 
fc,A 



'fe,A 



\ * 



g 



23 



E' 



23 



ttQED 



CtQED 



(,12 

"atom 



) 



p23 

"atom 



) 



(B.2) 



(B.3) 



(2vr)' Uk 

where oqed is the fine-structure constant. The transition matrix elements for the direct 
two-photon absorption (or emission) process, which is possible via the quadratic term 
OC A (ro + Sr,t), are zero when using the dipole approximation, as there would be no 
(momentum) operator left working on the electron states. Therefore it is necessary to 
look at higher orders of e'*"'^^ 

(ls| e*-^^ |3s) = (ls| [l + ik-dr-{k- 5rf /2 + O {{k ■ 5rf)] |3s) 



(ls| (fc • (5f)^ |3s) + O ((fc ■ (5f)^) 



(B.4) 



As for the hydrogen atom, we assumed (ls| dr |3s) = for parity reasons (Is and 3s 
are both spherically symmetric). Therefore we will take account of the quadratic order 
as the first non-vanishing order, while neglecting higher orders which should be much 
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smaller. This results in the following transition matrix element 
^i3 = il(i3|ei('^+^)-5^|3s) 



2m 

1 — (1-5 
4m 



(fc + fc) ■ 6r |3s) 

^ ^^,fc + fc)2(43„J% (B.5) 
and the following coupling constant (note that we have quadratic order in as well) 

4 (ZTT) m y/^k^k ^ ' 

For (yf^^ ^- , the atomic part of the coupling constant simply reads g^^ = / (2m) 

because there is no atomic transition occurring in the case of direct oc A {r^ + dr,t)- 
scattering. There is an additional factor 2 because the process arises as a mixed term 

^ 2 

in A (ro + Sr,t). Thus 

9ll .1 = 2 (gt* ■91) — = ^ ^— (ei ■ ei) . (B.7) 

Using the recently derived coupling constants, the Hamiltonian (HOj) can be 
rewritten in terms of atomic transition operators like = |ls) (2p| and 

photonic annihilation/creation operators (^j,J(^\)^ order to prepare the subsequent 
perturbation theory calculations. There, ( H3l) . we did not quote the full result, which 
reads 

k,\ 

/ ^ ^kk,\\ k,\ fc,A 



fc,fc,A.A 



'fefc,AA fe,A fe,A 

fe,fc,A,A 

fc,A ' 

+ y (7ii,,e-^(-^+^^+^»)*a, a--aiV(^+^)-^«+H.c. 

/ J ^kk,X\ fc,A fe,A 

fc,fc,A,A 

/ ^ ^fcfc,AA„ll fc,A fe,A 

fc,fc,A,A 

+ y J^e'^-^--i'+Els)t^f . ^^^^^^ 

/ ^ ^kk,X\ „11 fc,A fe,A 
fe,fc,A,A 

_ y ^MA^_i(^,+,-,)t i(fc+A,)..o + H.C.. (B.8) 
fe,fc,A.A 
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In section I3.H the remaining terms -ff{'j™(ro, t) were suppressed because they are 
unimportant for the subsequent analysis 

^iT(^o, t)= Yl 9lUe--^^-''^--^^^^^\ ^a, .^i3ei(^+^)-^o + h.c. 



fc,fc,A,A 



/ ^ ^kk,\Xgll k,X k,X 

k,k,\,X 



,13 



gll ^e^K-^+i^i3)tat a _ai3tgi(fc-fc).r.o 

^kk,\X „11 k,X k,X 



k,k,X,X ^ 
q^l - 

fc,fc,A,A 

^ 3 

In the perturbation theory calculations of the Hamiltonian (IB.Sp . all terms of order A 
or higher are ignored. Thus all the remaining terms (1B.9P can only contribute to the 
first-order amplitude. However, the first, third and fourth term all have a non-vanishing 
temporal phase (when a;i/2 < -E13) and thus vanish for infinite interaction time. The 
second term is not important as it starts from the 3s-level where in our initial state the 
atom occupies the Is-level. 

Appendix B.2. Integrating out the 2p-level 

The two-photon absorption probability of (IHHj) in standard first-order perturbation 
theory in case of S" = 1 is just P27 = T'^\9^ ^ where T is the total interaction 
time. Comparing this result with the two-photon absorption probability of the original 
Hamiltonian, f l55|) . one can extract the effective coupling constant g to 

I I ^ StT [g'fci^Ail kfc^.Aal 1 ^ 2aQED -^12-^23 ^atom m -j^qN 

T A A Tiy/uj^T A A ' y ■ ) 

Appendix C. Miscellaneous 
Appendix C.l. Commutation relations 

In order to simplify the results of the perturbation theory calculations, throughout the 
whole paper the following commutation relations (1C.4P and (1C.5|) were applied. They 
can be derived using the definition of a-^^^^y 

«l/2 = E (fc) V,^^^ , (C.l) 

k 

and the fundamental commutation relation for annihilation/creation operators 

6' (fc - fc) 6,-, . (C.2) 
From these two definitions, we can easily calculate the following commutation rule 



^kX^k'x 



\,A^V2 = ^V2^fc,A + /l/2(^)^AAi/2 " (^.3) 
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Applying (1C.3P repeatedly yields the commutation relations employed in section 13.21 



and 



fc,A 



+ /rm,A (4«,, + /2m.A. 



(C.4) 



(C.5) 



Appendix C.2. Scattering probability integration 

Here the derivation of the scattering probability integration is given in full detail. We 
start with expression (!57|) from section 13.41 where we already neglected the /^-terms 



fe,A,C=l,2 



E 



,12 I „12 Y 
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^kk,XXt 



^i(k—k)-ro 



(C.6) 



We then carry out the inner integration assuming constant integrand A; ^ fc^, which 
yields 



Pi, = 2(2vrf Yl 

fc,A,C=l,2 



11 



fcfeij,AA^ 



Ak, 



(C.7) 



Now we write out the outer integration over k explicitly in spherical coordinates such 
that 



k = k (sin 6 cos (p, sin 6 sin (p, cos ( 



bearing in mind that the integral over the radius is restricted to a small interval due to 
the Dirac delta function of the inner integration 



P,, = 2{2nf 



k(+Ak^ 



A,C=1,2 ki;~Akz 



k^dk / sin^d^ 
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The radius-integration is quickly resolved by again assuming a constant integrand, and 
what remains to be calculated is (remember that e^^ ■ -^^tom — ^atom) 



/»7r /»27r 

Ply = A{2TTf / sinede 
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QED^atom 
72 
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' ^12 %fcL,AA 
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- ^12 -"-atom 
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(C.IO) 



The two orthogonal, transverse polarization vectors can be written as 

ej[. = (sin^^, — COS0, 0) , 

el = (cos 6 cos 6, cos 6 sin — sin 9) . 



iC.U] 



Using the following integral equation which can easily be calculated and is valid for any 
vector V 

r2-K 



/"TT 1*2. 

V / sinOdO I 
^ Jo Jo 



V ■ e. 



xA|2 



-TT \V\ 



we arrive at the result already stated in section 13.41 (remember £. 
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